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ABSTRACT 

We present an eclipse timing variations analysis of the post-common envelope 
binary NSVS 14256825, which is composed by a sdOB star and a dM star in a 
close orbit (Porb ~ 0.110374 days). High-speed photometry of this system was 
performed between July, 2010 and August, 2012. Ten new mid-echpse times of 
NSVS 14256825 were analyzed together with all available eclipse times in the 
literature. We revisited the 0-C diagram using a hnear ephemeris and detected 
a clear orbital period variation. We investigated this variation on the assumption 
that it is a light travel time (LTT) effect. A full orbital motion analysis indicates 
that two LTT shifts with semi-amplitudes of ~ 20 s and ~ 5 s are superimposed 
on the 0-C diagram and are consistent with two circumbinary bodies. The 
best solution provides the orbital periods, Pc = 3.49 ± 0.21 years and P^ — 
6.86 ±0.25 years, and the projected semi-major axes, ac sin 7c = 1.9 ±0.3 AU and 
Qd sin Id — 2.9 ± 0.6 AU, for the circumbinary bodies. The masses of the external 
bodies are ~ 2.8 Mj^piter and ~ 8.1 Mj^piter, if we assume their orbits 
are coplanar with the close binary. Therefore NSVS 14256825 is composed by 
the close binary and two circumbinary planets in a 2:1 mean motion resonance. 
The closer planet in NSVS 14256825 has the minimum binary-planet separation 
among all known circumbinary planets in post-common envelope systems. 

Subject headings: planetary systems - binaries: eclipsing - binaries: close - 
subdwarf - stars: individual: NSVS 14256825. 



Based on observations carried out at the Observatorio do Pico dos Dias (OPD/LNA) in Brazil. 
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INTRODUCTION 



Planetary formation and evolution around binary systems have become important top- 



ics since the discovery of the first exoplanet around the binary pulsar PSRB1620-26 (Backer 



et al. 1993). Theoretical studies have indicated that circumbinary planets can be formed 
and survive for a long time (Moriwaki &; Nakagawa | |2004[ IQuintana fc Lissauer||2006D . Char- 
acterization of such planets in different evolutionary stages of the host binary is crucial to 
constrain and test the formation and evolution models. 

The common envelope (CE) phase in binary systems is dramatic for the planets survival. 



For a single star, Villaver & Livio (2007) have pointed out that planets more massive than 
two Jupiter masses around a main sequence star of 1 Mq survive the planetary nebula stage 
down to orbital distances of 3 AU. The CE phase in binary systems is more complex than the 
nebular stage in single stars and its interaction with planets is still an open topic. Besides, 
a second generation of planets can be formed from a disk originated by the ejected envelope 



(van Winckel et al. 2009 Perets 2010). Investigation of circumbinary planets in post-CE 
phase systems is fundamental to constrain observationally the minimum host binary-planet 
separation and to distinguish between planetary formation before and after the CE phase. 

To date, circumbinary planets have been discovered in seven eclipsing post-CE binaries. 
All those discoveries were done using the eclipse timing variation technique. The main 
features of these planets are summarized in Table |4j This table also shows the results on 
NSVS 14256825 presented in this paper. 

NSVS 14256825 (hereafter referred to as NSVS 1425) is an eclipsing post-CE binary and 



consists of a sdOB star plus a dM star with an orbital period of 0.110374 days (Almeida et 



al. 


2012 


1. It was discovered 


Survey i 


Wozniak et al. 


2004) 



NSVS 1425 is increasing at a rate of ~ 1.1 x 10~ s s . Beuermann et al. (2012b) presented 
additional eclipse times of NSVS 1425, suggesting from the analysis of the 0-C diagram the 
presence of a giant planet of ~ 12 Mjupitcr- 

In this study we present 10 new mid-eclipse times of NSVS 1425 obtained between July, 
2010 and August, 2012. We combined these data with previous measurements from the 
literature and performed a new orbital period variation analysis. In Section |2] we describe 
our data as well as the reduction procedure. The methodology used to obtain the eclipse 
times and the procedure to examine the orbital period variation are presented in Section [3j 
In Section H] we discuss our results. 
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2. OBSERVATIONS AND DATA REDUCTION 

The observations of NSVS 1425 are part of a program to search for echpse timing vari- 
ations in compact binaries. This project is being carried out using the facihties of the 
Observatdrio do Pico dos Dias, which is operated by the Laboratorio Nacional de Astrofisica 
in BraziL Photometric observations were performed using CCD cameras attached to the 
0.6-m and 1.6-m telescopes. Typically 100 bias frames and 30 dome flat-field images were 
collected each night to correct systematic effects from the CCD data. The photometric data 
are summarized in Table |2l 

The data reduction was performed using iRApj^tasks and consists of subtracting a master 
median bias image from each program image, and dividing the result by a normalized flat- 
field frame. Differential photometry was used to obtain the flux ratio between the target and 
a field star with constant flux. As the NSVS 1425 field is not crowded, the fluxes extraction 
were performed using aperture photometry. This procedure was repeated several times using 
different apertures and sky ring sizes to select the combination that provides the best signal- 
to-noise ratio. Figure [T] shows three normalized light curves folded on the NSVS 1425 orbital 
period. 



3. ANALYSIS AND RESULTS 



3.1. Eclipse fitting 



To obtain the mid-eclipse times for NSVS 1425, we generated model light curves using 
the Wilson-Devinney code (WDC) ( Wilson fc Devinney||1971 ) and searched for the best fit to 
the observed data. We used mode 2 of the WDC, which is appropriate to detached systems. 
The luminosity of each component was computed assuming stellar atmosphere radiation. 
The linear limb darkening coefficients, Xj, were used for both components. For the unfiltered 
light curves, V^-band limb darkening was used. The ranges of the geometrical and physical 



parameters (e.g., inclination, radii, temperatures and masses) obtained by Almeida et al. 



(2012) for NSVS 1425 were adopted as the search intervals in the fit. 



A method similar to the described in Almeida et al. (2012) was used for the fitting 



procedure. The WDC was used as a "function" to be optimized by the genetic algorithm 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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PIKAIA ( Charbonneau|l995 ). To measure the goodness of fit, we use the reduced Xred defined 
as 

2 1 I Oi — Ci 

i=i ^ * 

where Oi are the observed points, Ci are the corresponding models, ai is the uncertainty 
at each point, and n is the number of points. Figure [T] shows three eclipses of NSVS 1425 
superimposed with the best solution. To establish realistic uncertainties, we used the solution 
obtained by pikaia as input to a MCMC procedure (Gilks et al. 1996) and examined the 
marginal posterior distribution of probability of the parameters. The mid-eclipse times and 
corresponding uncertainties were obtained from the median value of the marginal distribution 
of the fitted times and the 1-a uncertainties from the corresponding 68% area under the 
distribution. The results are presented in Table |4} 



3.2. Linear ephemeris 

To determine an ephemeris for the NSVS 1425 orbital period, we analyzed our mea- 
surements together with all available eclipse times in the literature after converting them to 
barycentric dynamical time (TDB). Table |4] shows all eclipse times available for NSVS 1425. 
Fitting the data using a linear ephemeris, T^nm = Tq + E x P^in, we obtain, 

Trnin = TDB 2454274.2086(1) + 0.110374165(1) x E, (2) 

where T^in are the predicted eclipse times, Tq is a fiducial epoch, E is the cycle count from 
To and Pbin is the binary orbital period. The best fit yields a Xred ~ ^6. The residuals of the 
observed times with respect to Eq. [2] are shown in the 0-C diagram of Fig. |2] 



3.3. Eclipse timing variation 

Fig. |2] shows that a linear ephemeris is far from predicting correctly the NSVS 1425 
echpse times. The large value of Xred suggests the presence of additional signals in the 0-C 
diagram, which can be explained by the light travel time (LTT) effect. 

The LTT effect will show up as a periodic variation in the observed eclipse times when 
the distance from the binary to the observer varies due to gravitational interaction among 



the inner binary and an external body (Irwin 1952). To fit the NSVS 1425 eclipse times 



taking this effect into account, we used the following equation, 

n 

Tmin = To + E X Pun + J2 (3) 



n=l 
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where 

(4) 

is the LTT effect. In the last equation, Ki = ai sin Jj/c is the time semi-amphtude, Cj is the 
eccentricity, Ui is the argument of periastron, and fi is the true anomaly. These parameters 
are relative to the orbit of the inner binary center of mass around the common center of 
mass consisting of the inner binary and of the i—th planet. The parameters Oj, and c in 
the semi-amplitude equation are the semi-major axis, the inclination, and the speed of light, 
respectively. Notice that we do not consider mutual interaction between external bodies in 
this analysis. 

Initially we fitted Eq. [sjto the data with only one LTT effect. The resulting Xred dropped 
to 6.8, but the new residuals showed evidences of another cyclic variation. Adding one more 
LTT effect in Eq. |3| the resulting x^ed improves to L85. The pikaia algorithm was used 
to search for the global optimal solution, followed by a MCMC procedure to sample the 
parameters of Eq. |3] around the best solution. Figure |3] shows the resulting 0-C diagram 
and Table |4] shows the best fit parameters with the associated ±68% uncertainties. 



Zi, 

c 



l-e 



1 + Ci cos fi 



sin(/i + Ui 



4. DISCUSSION AND CONCLUSION 

We have analyzed 10 new mid-eclipse times of the post-CE binary NSVS 1425. The 
0-C diagram shows a complex orbital period variation. We investigate the properties of this 
variation under the assumption that it is caused by the LTT effect. Our analysis indicates 
that two distinct effects can be traced out in the 0-C diagram, one with amplitude ~ 20 s and 
the other with amplitude 5 s. The associated orbital periods correspond to 6.86 ± 0.25 years 
and 3.49 ± 0.21 years, respectively. Despite the partial coverage (~ 75%) of the longer 
cycle, the indicative resonant motion with a period ratio of 2:1 gives more confidence to the 
circumbinary bodies hypothesis. 

Assuming two external bodies around NSVS 1425 and using the close binary mass Mbin = 



O.528M0 (Almeida et al. 2012) the lower mass limit for the two circumbinary bodies are 
Mc sin Jc ~ 2.8 Mjupiter (inner body) and M^sin/d ~ 8.0 Mjupiter (external body). The close 
binary has an orbital inclination of 82.5 degrees ( Almeida et al.|[2012 ). Under the additional 



assumption of coplanarity between the two external bodies and the inner binary, NSVS 1425 c 
and NSVS 1425 d would both be giant planets with Mc ~ 2.8 Mjupiter and ~ 8.1 Mjupiter- 

NSVS 1425 is the eighth post-CE system with circumbinary planets and the fourth 
system with two planets (see Table |4]). In such systems, there are two principal scenarios for 
planetary formation: (i) first generation planets formed in a circumbinary protoplanetary 
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disk; and (ii) second generation planets originated from a disk formed by the ejected envelope 



(van Winckel et al. 2009 Perets 2010) 



Considering the first scenario, could the two circumbinary planets in NSVS 1425 survive 



the CE phase? Bear & Soker (2011) estimated the orbital separation between the progenitor 



of an extreme horizontal branch (EHB) star (sdB or sdOB) and a planet before the CE phase 
by the equation, gq ~ Mehb^ehb/^pto, where Mehb is the mass of the sdOB star, and oehb 
is the present orbital separation of the planet. Assuming that the progenitor of the sdOB 
star in NSVS 1425 had a mass of Mpj-o = 1.0 Mq, and neglecting accretion by the companion, 
the orbital separation of the planets before the CE would be aoc — 0.8 AU and ood — 1.3 



AU. For a single star, Han et al. (2002) pointed out that the maximum radius at the tip of 
the EHB (-Rehb) is ~ 0.8 AU. Hence, the inner planet is on the verge of being engulfed by 
the CE. Moreover, the tidal interaction causes a planet to spiral inwards if its orbital radius 



is smaller than < 3 -Rehb (Villaver & Livio 2007). Therefore the two circumbinary bodies 



in NSVS 1425 would not survive the CE phase. On the other hand, Taam & Ricker (2010) 



showed that the CE size can be much reduced if the EHB star is part of a binary system, 
because once the secondary is engulfed by the envelope, the CE is totally ejected in only 
~ 10^ days, stopping the envelope expansion. Therefore, the maximum radius of the CE is 
around the initial distance between the close binary components. Thus, if the close binary 
separation in NSVS 1425 before the CE phase was < 0.27 AU, the two planets could survive 
the CE phase. 

For the second scenario, the principal question to investigate is: was there time enough 



after the CE phase to form giant planets? Kley ( 1999 ) showed that the typical time scale to 
form giant planets in protostellar disks is ~ 10^ years. The lifetime of a binary in the EHB 
phase is ~ 10*^ years (Dorman et al. 1993 Heber 2009). As the NSVS 1425 primary star is 



in the post-EHB phase (Almeida et al. 2012), we conclude that there was time enough to 



form the two circumbinary planets after the CE phase. Therefore the second generation of 
planets is also a viable scenario for NSVS 1425. 

Finally, among all known candidates to be circumbinary planets in post-CE systems, 
the inner planet in NSVS 1425 has the minimum binary-planet separation, OcSin = 1.9 ± 
0.3 AU. 
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Fig. 1. — Primary eclipse of NSVS 1425 folded on the binary orbital period. The solid line 
represents the best fitting performed with the Wilson-Devinney code (see the Section [s]). 
The upper and lower light curves were displaced vertically 0.2 units for better visualization. 
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Fig. 2. — O— C diagram of the eclipse timings of NSVS 1425 built using Eq. |2} Our data are 
presented with full squares. 
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Fig. 3. — The upper panel shows the (0-C) diagram of the echpse times of NSVS 1425 built 
with respect to the linear part of the ephemeris in Eq. [3j Our data are presented with full 
squares and the solid line represents the best fit including the two LTT effects. The second 
and third panels display separately the two LTT effects [ti and T2). The lower panel shows 
the residuals around the combined fit. 
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Table 1: Circumbinary planets discovered in post-common envelope systems. 



Name 



orb 



Mp sin / 



a sin/ 



Spec. Type References 





days 


A/f, 

Jupiter 


ATT 








IN O V O L^Zoyi\D JC 


LZ 1 


9 8 
Z.o 


1 Q 

i.y 


n n 
u.u 


crlOR-LrlA/T 
SQ^JlJ-rQiVi 


+ r^lO 0"riT /H T T 

xnis btU-uy 


U Zj rOTyPi-D )C 


1 ni 7 
iyi ( 


7 7 
( . ( 


9 fi 

z.o 


U.UO 


UA+aivi (^r olar ) 


i 


HU Aqr(AB)c 


2226 


4.5 


3.32 


0.11 


DA+dM (Polar) 


2 


NSVS 1425(AB)d 


2506 


8.0 


2.9 


0.52 


sdOB+dM 


this study 


NN Ser(AB)c 


2605 


4.0 


3.2 


0.05 


DA+dM 


3 


NY Vir(AB)c 


2900 


2.3 


3.3 




sdB+dM 


4 


RR Cae(AB)c 


4346 


4.2 


5.3 


0.0 


DA+dM 


5 


HW Vir(AB)c 


4640 


14.0 


4.69 


0.4 


sdB+dM 


6 


HU Aqr(AB)d 


5155 


5.7 


5.81 


0.04 


DA+dM (Polar) 


2 


NN Ser(AB)d 


5571 


6.71 


5.32 


0.22 


DA+dM 


3 


UZ For(AB)d 


5844 


6.3 


5.9 


0.04 


DA+dM (Polar) 


1 


DP Leo(AB)c 


10227 


6.05 


8.18 


0.39 


DA+dM (Polar) 


7 



References. — (1) Potter et al. (2011); (2) Hinse et al. (2012); (3) Horner et al. (2012); (4) Qianetal 



(2012a); (5) Qian et al. (2012b); (6) Beuermann et al. (2012a|; (7) Beuermann et al. (2011) 



Table 2: Log of the photometric observations 



Date 




^exp(s) 


Telescope 


Filter 


2010 Jul 30 


300 


20 


0.6-m 


Rc 


2010 Jul 31 


450 


20 


0.6-m 


Rc 


2010 Aug 18 


800 


10 


0.6-m 


Rc 


2010 Nov 20 


350 


20 


0.6-m 


Ic 


2011 Jul 06 


1255 


2 


1.6-m 


Ic 


2011 Jul 07 


1300 


1 


1.6-m 


Clear 


2011 Aug 06 


435 


5 


1.6-m 


V 


2012 Apr 24 


1550 


1.5 


0.6-m 


Clear 


2012 May 30 


600 


3 


0.6-m 


Clear 


2012 Aug 12 


1330 


2 


1.6-m 
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Table 3. Eclipse times for NSVS 1425 



Cycle 


Time (BJD TDB) 


C (s) Eclipse 


Ref. 


± 


Z^:O^Z 1 4.ZUoo(^i ) 


Q n T 
y.u 1 


1 
1 


70 


0/|c;/l080 '\KKQ(0\ 
Z404ZoZ.i00y i^Z ) 


1 Q /I T 

iy.4 1 


i 


/ o 


0/1 1^/1080 0RR1 fO\ 
Z404ZoZ.ZDDi yZ ) 


/I Q T 
4.0 1 


i 


iUo 


0/1 1^/I08fi 1 0Q1 ('1 ^ 
Z404ZoD.iZyH^i J 


/in T 

-4.U 1 


± 


1 70 


O/ICi/IOQ*? 1QQ0('1^ 

Z4o4zyo.iyoz^^± ) 


Q T 

y.z 1 


i 


ioU 


Z404Zy4.U ( uzyL ) 


Q 8 T 

y.o i 


1 
i 


ioi 


Z4D4Zy4.ioDD(^i ) 


10 1 T 
iZ. i i 


1 
i 


1 QD 

iyu 


z^oHiZvo.i 1 yy i^i ) 


R T 
D.Z i 


1 
i 


oiu 


Z404oUy.Uo ) 


Q T 
Z.O i 


1 
i 


Q1 7 


Z404ouy.iy ( oyi ) 


/II T 
-4. 1 i 


1 
i 


QOc; 
OZO 


OA c;/iQi n nan/i \ 

Z404oiU.UoU4|^i ) 


1^ 1 T 
O.i 1 


i 


QKO 
ODZ 


OA ^A'^T A 1 RAO(^ \ 
Z404oi4.iD4Z^^± ) 


1 Q T 
i.O 1 


i 


OoU 


Z404OlD.10Uyi 1 ) 


17 T 
-L.I i 


1 


QQ7 


0/ii^/iQi8 ^\07A^^\ 

Z404oiO.UZ ( 4(^iJ 


1 n Q T 
iU.O 1 


± 


4Uu 


0/ic;/iQiQ nonfi^'i^ 

Z4D4oiy.UZUD\^i j 


/I T 
-4.Z i 


i 


Ar\7 


OA^A'i^Q 1 'il o^'i ^ 

Z4040 iy . id iZ i j 


^ ^ A T 


1 
i 


/I/I Q 


OA c;/iQOQ 1 n/1 \ 

Z404oZo.iU40(^i ) 


n 7 T 


1 
i 


/I KO 
40Z 


OAKA'^OA ^\Q7Q(^^ 

Z404oz4.uy ( y i 1 


T 

O.O i 


1 


fiQO 
ooz 


Z404oDD.U4Ui(^i ) 


A Ck T 

4.y 1 


1 
i 


lUlO 


0/l'^/l'^8R ^f{Q'^(«\ 
Z4040oD.ODyo^^O ) 


00 fi T 
-zy.D i 



z 




Z404DoD.D / Oy ) 


8 8 T 

-O.O i 



z 




Z^UUUO / .00U04i Z J 


■^4 n T 

-04.U 1 



z 


7037 


2455050.91137(2) 


-33.7 I 


2 


7304 


2455080.38128(7) 


-33.0 I 


2 


7322 


2455082.36800(2) 


-34.1 II 


2 


9823.5 


2455358.46897(5) 


-35.0 I 


2 


9959 


2455373.42474(7) 


-29.1 I 


2 


10131 


2455392.40910(2) 


-29.0 I 


2 


10279 


2455408.74442(2) 


-33.7 I 


this study 


10287 


2455409.62744(2) 


-31.7 I 


this study 


10451 


2455427.72877(1) 


-34.0 I 


this study 
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Table 3 — Continued 



Cycle 


Time (BJD-TDB) 


0-C (s) 


Eclipse 


Ref. 




245544Q 251 76('5'l 


-31 Q 




'X 




infi47 


9455440 '?6915f9'l 


-30 6 




'X 






9455459 9'?18Ql'9'l 

ZjrtOOrtO^.^O J-(Jiy\^ I 


-99 6 






1 1 1 4fi 


9455504 4Q405l'1'l 


-30 4 


TT 

XX 


"I" n 1 G G"!" n H AT" 


1976'? 


9455689 Q1400f9'l 


-91 9 







197QQ 


9455686 88745l'9'l 


-99 7 







19700 c; 


9455686 Q4970('S'l 


-17 4 


TT 

J. J. 


9 




945571 7 57146l'l 


-93 5 

ZjO. 




"X 




j.OO\JO 


9455740 6Q0S61 ((\\ 


-99 




"1" ill G G"!"!! f\ \T 
tlilo OLLlLlj 


1 3377 


9455750 68371 7('4'l 


-93 






1 3fi9Q 


9455778 4Q8061 


-18 3 

xo.o 






1 34fiQ 


9455760 83818(''?'l 


-19 6 




9 


1 34fiQ 5 


9455760 8Q340('3'l 


-17 1 

X 1 . X 


TT 

XX 


9 


1 '^470 


9455760 Q4855f1 1 


-90 

Z/U. u 




9 




9455769 Q'?5S9l'9'l 
jLi'-too i xjiLi.ijooo^y^ J 


X 1 . 




9 


LOO X X 


9455765 47'?87('9'l 
jlj'^oo i \jo i 00 i \^ ) 


-99 




9 


1 S^49 


9455768 8Q54Q('4'l 


-90 

ZjU. U 




9 


1 3639 


2455778 82Q1 5^1 1 


-90 q 




9 


1 3689 


9455784 34781 {1\ 


-95 9 




9 


1 3768 
xo / uo 


94557Q3 84006n 


-19 9 




9 


1 3897 


9455800 35917('9'l 


-16 5 




9 


1 3898 


2455800 46251 (2) 


-19.2 




2 


13845 


2455802.33887(5) 


-19.4 




2 


13846 


2455802.44920(4) 


-23.4 




2 


13872 


2455805.31896(2) 


-20.5 




2 


13873 


2455805.42932(2) 


-21.5 




2 


13899 


2455808.29907(2) 


-20.2 




2 


14062 


2455826.29008(2) 


-18.0 




3 


14089 


2455829.27017(2) 


-19.1 




3 


14379 


2455861.27873(2) 


-14.6 




3 


14397 


2455863.26542(3) 


-18.5 




3 
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Table 3 — Continued 



Cycle 


Time (BJD-TDB) 


0-C (s) Eclipse 


Ref. 


14400 


2455863.59656(1) 


-17.1 I 


2 


14406 


2455864.25879(3) 


-18.3 I 


3 


16024 


2456042.844216(4) 


-16.4 I 


this study 


16350 


2456078.826240(8) 


-14.4 I 


this study 


17019 


2456152.666554(6) 


-12.6 I 


this study 



References. — (1) |Wils et aL (2007); (2) [Beuermann et alT 
( |2012a[ ); (3) |Kilkenny fc Koen| ( |20l2| . 
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Table 4: Parameters for the linear plus two-LTT ephemeris of NSVS 1425. 



Linear ephemeris 


Parameter 


Value 


Unit 


-Pbin 


0.1103741681(5) 


days 


To 


2454274.20874(4) 


BJD(TDB) 


Ti term 


Parameter 


Value 


Unit 


P 


6.86 ±0.25 


years 


T 


2456643 ± 60 


BJD(TDB) 


Obin sin / 


0.042 ± 0.004 


AU 


e 


0.52 ± 0.06 






97.5 ±8 


degrees 


f{M) 


(1.6 ±0.5) X 10-6 


Mo 


Msin/ 


8.0 ±0.8 


-^Jupiter 


a sin/ 


2.9 ±0.6 


AU 


T2 term 


Parameter 


Value 


Unit 


P 


3.49 ±0.21 


years 


T 


2455515 ± 55 


BJD(TDB) 


Obin sin I 


0.0099 ± 0.0006 


AU 


e 


0.00 ±0.08 




CO 


11.4 ± 7 


degrees 


f{M) 


(8.0 ±2.0) X 10"^ 


Mo 


Msin/ 


2.8 ±0.3 


-^Jupiter 


a sin/ 


1.9 ±0.3 


AU 


Xred 


1.85 





